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A new method to compute observables at many values of the parameters A for a model with lattice action 
S((f>, A) is described. After fixing a reference set A r of parameters, a single simulation is carried out by using a 
"reference action" S((j> r , A r ) to generate configurations of the field </> r . Then a suitable analytic transformation is 
performed from the configurations of (jf to the ones corresponding to the action S(<j>, A). Such a transformation 
allows to obtain the observables for values of the parameters A close to A r . I present studies on the reliability of 
the algorithm in the case of the cf> 4 model in 2 dimensions. 



1. INTRODUCTION 

Usually in Monte Carlo (MC) simulations ob- 
servables are measured at the values of the La- 
grangian parameters at which the run is carried 
out. However, it is often necessary to evaluate 
observables at many different values of the pa- 
rameters, for instance, in problems requiring a 
fine tuning. This poses dramatic limitations to 
the numerical studies. 

The reweighting technique Q has long been a 
way out of this problem. Let us consider a single 
MC run performed according to a reference action 
S[4>, \ r ] with given values X r for the Lagrangian 
parameters. Field configurations are generated 
with probability e~ s ^' X ^ and expectation values 
for observables are computed averaging over these 
configurations. Since the equilibrium distribution 
is known, physical quantities can be estimated for 
different parameters without performing a new 
simulation. The expectation value of an observ- 
able for the set of parameters A can be calculated 
as a reweighted average from the configurations 
generated by the MC simulation performed at the 
set A r . Moreover, the field configurations gener- 
ated according to e~ 5 ^ ,Ar l must be a good sam- 
pling for the distribution e~ s ^'^. In the limit 
of an infinite number of configurations, e~ s ^'^ 
is always correctly reconstructed but, as this is 
not the real case, there is a systematic distortion 
of the equilibrium distribution. To minimize this 
effect, the parameter set A should not differ much 
from the reference set A r . 



Recently, an alternative technique B has been 
proposed to compute observables at various val- 
ues of the couplings and masses by a sin- 
gle MC simulation. This method is based on 
a "field transformation" from a reference field 
<fi r , distributed with probability proportional to 
e -s[<p r ,x r ]^ £ Q a -ggi^ ^ distributed according to 
e -5[0,A] rp ne fc rans f orma tion is defined by the 
equation 



Vcj) e 



(1) 



Thus, an importance sampling for the field <jf 
can be transformed into an importance sampling 
for the field <j>. I will explore this technique by 
using a simplified version of ([!]). More precisely, 
I consider only the terms of the action depending 
on single sites. This simplification implies that it 
is necessary to reweight the data for observables 
with a suitable remainder action 8S[(f) r : A, X r ]. 
Work is in progress to take into account the full 
action. 

Here I study how reliable is the field transfor- 
mation technique compared with the reweighting 
method. The model I consider is the lattice </> 4 
scalar field theory in two dimensions. 



2. THE CASE STUDIED: THE TWO- 
DIMENSIONAL </> 4 MODEL 

I consider the 2 — d lattice scalar <f> 4 field theory 
with periodic boundary conditions and L 2 sites. 
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The lattice action is 
«S[<M = -^</»„(0„+p+0„- Al )+^« A (0„)(2) 



where ( A = (A2, A4) ) 
v\<t> n ) = )-(\ 2 +A)4>l + )\ 



(3) 



The sum over fi is to be intended in the positive 
directions. 

The expectation value for a monomial P[cj>] = 
in the field is defined to be 



-"ni Vn 2 ■ ■ ■ Yrik 



< P[4>] > x ^ I / [J dK e- s ^P[4>] (4) 

n 

where Z is the partition function. 

A generic field transformation <f> — 4>{<t> r ) is de- 
fined by the jacobian matrix: 



Unra — 



(5) 



Re-expressing the field polynomial P[4>] in terms 
of 4> r , we obtain a function V\(jf \ of the reference 
field. The same transformation in the measure of 
the functional integral yields 

JJd0 B e- 5 ^ A I = JJdft e -^-M+ssi^.xxn {6) 



where SS, the "remainder action", is 

SS[cj) r , A, A r ] = -S[<t>, A] + S[cb r , A r ] + Tr In J (7) 

Now the Green function assumes the following 
form 



< PW\ >> 



< V[(j) r } e 5S ^ > X r 
< e ss[4> r ] >Xr 



(8) 



The remainder action vanishes if (Q) is exactly 
solved, otherwise it enters as a reweighting term. 
The integration of (||) for a pair A = (A2,A4) al- 
lows to compute V[<f> r ] and 5 t S[0 r ,A, A r ] for each 
MC generated configuration of the reference field 
(j) r . Then by averaging over the reference field 
configurations, one obtains the Green function 
< P[<p] >\ corresponding to the Lagrangian with 
parameters A. 



2.1. Form of the Field Transformation 

The field transformation defined by (|l|) is not 
easily integrable. Thus, a simpler case has been 
considered. In the simulation, I have used the 
following diagonal jacobian: 



J n 



C5 n 



(9) 



where C is a constant. The solution of the 
previous differential equation with the condition 
<b n — when 6Z. = is 



dcjy e~ v w = C 
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-v x (0) 



(10) 



The fields (f> and <p r are non-compact and the con- 
stant C is determined requiring that tfi n — > 00 
when (j)^ — > 00. 

With the choice (||), it follows 



Tr In J = 5> A (0„)-^(<O 



(11) 



and consequently a non vanishing remainder ac- 
tion has to be used in (ph to evaluate Green func- 
tions as a ratio of reweighted averages. 

2.2. Field Transformation and Reweight- 
ing 

In this section I report on the tests of con- 
fidence performed for the field transformation 
method. These results are compared with those 
obtained by the reweighting technique. This anal- 
ysis extends the first studies of reliability carried 
out for the (j> 4 scalar field theory in 4 dimensions 
I have examined the susceptibility and the 
expectation value of the action S[<f>, A]. 

The MC simulation has been performed at the 
reference point A r = (—.2, .23) and 66 pairs of 
parameters have been considered: 28 pairs had 
A2 = A2 (group 1), while the remaining 38 ones 
had A4 = A4 (group 2). Moreover, independent 
MC simulations have been carried out to check 
the results. 

In the reference simulation 210000 uncorrelated 
data for the action and 40000 uncorrelated data 
have been collected for the susceptibility; the MC 
test runs have amounted to 40000 uncorrelated 
data. 
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The figures below show the results for the ex- 
pectation value of the action. 
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Figure 1. Expectation value of the action for the 
group 1 pairs 



1940.0 



1920.0 



1900.0 



1880.0 



1860.0 



1840.0 



-3^ 



field transformation 

* reweighting 
Preference simulation 

• MC test runs 



-0.30 -0.26 -0.22 -0.16 



-0.14 -0.10 



Figure 2. Expectation value of the action for the 
group 2 pairs 

Both for the group 1 and the group 2 pairs, 
the reweighting technique, the field transforma- 
tion method and the MC test runs agree within 
the error bars, except at small values of A2 and A4 
where the field transformation method performs 
better. 

The results for the susceptibility are analogous: 
the agreement between reweighting, field trans- 
formation and MC test runs is within one stan- 
dard deviation provided that the values of A2 and 
A4 are not small and the critical region is not ap- 
proached. 



As a consequence of the choice (j^) for the 
field transformation, the statistical error analysis 
is similar to that performed for the reweighting 
technique ||. 

3. CONCLUSIONS 

The field transformation method is a general- 
ization of the reweighting technique. The simple 
case considered of a field transformation depend- 
ing on single sites does not provide a substantial 
improvement respect to reweighting. Though the 
studies presented are preliminary, the results ob- 
tained are similar to those given by the reweight- 
ing technique. However, including the full action 
in the field transformation could yield a consis- 
tent improvement in the algorithm. Work is in 
progress in this direction. The aim is to modify 
the field transformation in order to eliminate the 
need of the weight. 
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